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The best level of ordering and straightening of carbon nanotube arrays is often achieved when they 
are grown in a dielectric matrix, so such structures present the most suitable candidates for future 
channeling experiments with carbon nanotubes. Consequently, we investigate here how the dynamic 
polarization of carbon valence electrons in the presence of various surrounding dielectric media 
affects the angular distributions of protons channeled through (11, 9) single- wall carbon nanotubes. 
Proton speeds between 3 and 10 a.u., corresponding to energies of 0.223 and 2.49 MeV, are chosen 
with the nanotube's length varied between 0.1 and 1 /^m. We describe the repulsive interaction 
between a proton and the nanotube's atoms in a continuum-potential approximation based on the 
Doyle- Turner potential, whereas the attractive image force on a proton is calculated using a two- 
dimensional hydrodynamic model for the dynamic response of the nanotube valence electrons, while 
assigning to the surrounding medium an appropriate (frequency dependent) dielectric function. The 
angular distributions of channeled protons are generated using a computer simulation method which 
solves the proton equations of motion in the transverse plane numerically. Our analysis shows that 
the presence of a dielectric medium can strongly affect both the appearance and positions of maxima 
in the angular distributions of channeled protons. 

PACS numbers: 61.85.-|-p, 41.75.Ht, 61.82.Rx, 79.20.Rf 
Keywords: nanotubes, channeling, dynamic polarization, rainbows. 



I. INTRODUCTION 



While progress in theoretical modeling and computer 
simulations of ion channeling through carbon nanotubes 
has reached a mature level, as recently reviewed in Ref. 
[l], efforts of experimentalists have only recently begun 
to bear fruit in this important research area. Because 
the best level of ordering and straightening of carbon 
nanotubes is achieved when they are grown in a dielec- 
tric matrix, such structures are perhaps the most suitable 
candidates for ion channeling through carbon nanotubes. 
It thus came as no surprise when Zhu et al. recently 
reported the first experimental data on He"'" channeling 
through an array of well ordered, multi-wall carbon nan- 
otubes (MWNTs) which were grown in a porous anodic 
aluminum oxide (AI2O3) membrane. On the other hand, 
carbon nanotubes have also been grown selectively within 
etched ion tracks in Si02 layers on Si by another exper- 
imental group 3], thus offering an interesting possibil- 
ity for realization of ion channeling through individual, 
single-wall carbon nanotubes (SWNTs) at a wide range 
of ion energies. In addition, in many applications of car- 
bon nanotubes it is desirable to have them embedded in 
a dielectric such as Si02 @|, or clamped by a metal shield 
[l| made of nickel, which is known to most readily bind 
to carbon nanotubes Q. 
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For ion channeling at the low (keV) and high (GeV) 
ends of the energy range, the surrounding material would 
predominantly serve as a passive container of carbon nan- 
otubes when the dynamics of ion motion is concerned. 
However, ions moving with medium (MeV) energies will 
induce strong dynamic polarization of valence electrons 
in the nanotubes which in turn will give rise to a size- 
able image force on the ions, as well as a considerable 
energy loss due to the collective, or plasma, electron ex- 
citations. The dynamic image force has been recently 
shown to give rise to the rainbow effect in the angular 
distributions of protons channeled through short (11, 9) 
single-wall Q and double-wall carbon nanotubes in 
free space, which is not otherwise observable in simula- 
tions of ion channeling through chiral carbon nanotubes 
using the continuum approximation for the interaction 
potential [ll|. Obviously, the presence of dielectric me- 
dia may affect these dynamic polarization forces, as well 
as the resulting ion trajectories, making the analysis of 
such effects in ion channeling through carbon nanotubes 
a timely task, which we take up in this contribution. In 
particular, we are interested here primarily in the effects 
of dielectric media on the dynamic image force on chan- 
neled ions and their angular distributions. 

The rainbow effect occurs and plays an important role 
in photon scattering from water droplets , nucleus- 
nucleus collisions flO. [Ill , [l^ , atom or ion collisions with 
atoms or molecules [ij] , electron- molecule collisions 1J| , 
atom or electron scattering from crystal surfaces [islll^. 
and ion channeling in crystals [itI, [l^. Moreover, the 
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rainbow effect lias been investigated recently in the con- 
text of grazing scattering of atoms from metal surfaces 
under channeling conditions by Schiiller et al. [l^ who 
showed that precise measurements of the well-defined 
maxima in the angular distributions of scattered atoms, 
attributed to the rainbow effect, can give detailed infor- 
mation on the interaction potential of the atoms with 
the metal surfaces. On the other hand, the theory of 
crystal rainbows has been formulated as the proper the- 
ory of ion channeling in thin crystals [l^, and has been 
subsequently applied to ion channeling in short carbon 
nanotubes [2l|, (22|, |2^, [23| • It is therefore expected that, 
in analogy with the surface channeling experiments , 
measurements of the rainbow effect in carbon nanotubes 
may give precise information on both the atomic config- 
uration and the interaction potentials within such struc- 
tures, which are not completely known at present. 

Previously reported simulations of ion channeling in 
carbon nanotubes [1^ [H, [23| paid virtually no attention 
to the effects of dynamic polarization of the nanotube 
valence electrons. However, this process is expected to 
contribute to the ion energy loss and to give rise to strong 
image forces on the medium-energy ions [2^, as was re- 
cently demonstrated in the computer simulations of an- 
gular distributions of protons channeled through chiral 
single- walled carbon nanotubes in vacuum (29j . The im- 
portance of the image force has also been emphasized 
in the related area of ion transmission through cylindri- 
cal channels in metals [13, HH, IH, [H, [s^l . Whereas in- 
teresting parallels can be drawn between ion channeling 
through carbon nanotubes and ion transmission through 
capillaries, it is important to notice several crucial differ- 
ences. Namely, while the lengths of capillaries in met- 
als can be comparable to those of nanotubes consid- 
ered here for ion channeling, their diameters are typi- 
cally an order of magnitude, or more, larger than those 
of carbon nanotubes, making their aspect ratios consid- 
erably smaller. More importantly, the inner surfaces of 
such broad channels in metals usually appear to be quite 
rough, which does not seem to affect too much the trans- 
mission of slow, highly-charged ions through such struc- 
tures 1121. However, it is questionable whether con- 
ditions required for channeling of fast ions leading to the 
rainbow effect can be met in capillaries with such char- 
acteristics [l^ H^. On the other hand, in cases where 
carbon nanotubes are grown in amorphous channels in a 
dielectric such as AI2O3 [g,] and Si02 or are coated 
by an amorphous layer of metal [1, , it is precisely the 
regular atomic structure of carbon nanotube that acts as 
a smooth " sleeve" , or " mantle" covering the underlying 
rough surface of the surrounding material, thus enabling 
ion channeling through such structures. It remains to be 
seen, however, whether large, multi-walled carbon nan- 
otubes can be grown inside the broad capillaries in met- 
als enabling some sort of ion channeling in their interior 
hollow regions. 

We shall investigate here how dynamic polarization of 
carbon valence electrons influences the angular distribu- 



tions of protons channeled in (11, 9) SWNTs with differ- 
ent dielectric media surrounding the nanotubes. We con- 
sider proton speeds between 3 and 10 a.u., and nanotube 
lengths between 0.1 and 1 //m. The image force acting 
on an ion moving in a nanotube surrounded by a dielec- 
tric medium has been recently calculated by means of a 
two dimensional (2D) hydrodynamic model for the car- 
bon valence electrons, while the surrounding medium was 
described by a suitable frequency dependent dielectric 
function [U, HI]. In the present simulations, we shall use 
the van der Waals radius of a carbon atom (0.17 nm) to 
approximate the distance between the nanotube wall and 
various dielectric media [1^, . This value also agrees 
with our density functional theory (DFT) calculations of 
the average equilibrium separation between graphene and 
a Ni (111) surface, based on the methodology described 
in Ref. 41. Other details of our simulation are similar to 
those reported earlier [3, 0| • 

After outlining the basic theory used in modeling the 
dynamic polarization effects of carbon nanotubes, we 
shall discuss the results of our ion trajectory simulations 
and give our concluding remarks. Atomic units will be 
used throughout unless explicitly stated otherwise. 

II. THEORY 

The system under investigation is a proton moving 
through an (11, 9) SWNT surrounded by a dielectric 
medium. The z-axis coincides with the nanotube axis 
and the origin lies in its entrance plane. The initial pro- 
ton velocity vector, w, is taken to be parallel to the z-axis. 
Proton speeds between 3 and 10 a.u., corresponding to 
energies of 0.223 and 2.49 MeV, are considered, with the 
nanotube's length varied between 0.1 and 1 /im. 

We also assume that the repulsive interaction between 
the proton and the nanotube atoms may be treated clas- 
sically, using the Doyle- Turner expression for the proton- 
carbon atom interaction potential [i^ . averaged axially 
and azimuthally [44]. The repulsive potential for 
proton channeling through the nanotube is then of the 
form 

167rdZiZ2 

4 

j=i 

where Zi = 1 and Z2 — Q are the atomic numbers of pro- 
ton and carbon atoms, respectively, d = 2a is the nan- 
otube diameter, I is the C-C bond length, r is the distance 
between the proton and nanotube axis, /q is the modified 
Bessel function, and aj = {0.115,0.188,0.072,0.020} and 
6j = {0.547, 0.989, 1.982, 5.656} are fitting parameters in 
atomic units [43 |. 

The dynamic polarization of the nanotube is treated 
by a 2D hydrodynamic model of the nanotube valence 
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electrons, based on a jellium-like description of the ion 
cores on the nanotube wall [H, [s^] ■ This model includes 
both axial and azimuthal averaging consistent with our 
treatment of the repulsive interaction. The self energy, 
or the image potential, i?soif , for a single ion of charge Z\ 
at position ro(^) is defined by 

-B^elf = (^l/2)*ind(ro(t),t), (2) 

where $ind is the potential at the point r , given in 
cylindrical coordinates by r = {r, z}, which is induced 
in the system by the presence of ion. After performing 
the Fourier transform with respect to time, and follow- 
ing the method of Doerr and Yu [i^, we consider the 
total electric potential to be the sum of the external per- 
turbing potential, ^cxtl?', w), and the induced potential, 
*i'ind(?', w), due to polarization of the nanotube and the 
dielectric boundary by the proton, so that 

4>(f,w) = $ext(^t^) + $ind(r,t^)- (3) 

The Poisson equation then gives 

V^$oxt(r, = -47rpext(r, w), (4) 



The Fourier transform of the external perturbing po- 
tential due to a single proton of charge Z\ = \ moving 
parallel to the nanotube axis with constant speed u, such 
that ro(<) = {ro,(^o,wi}, is given by 

where £nt is the background optical dielectric constant 
for the nanotube (for which we use Ent = 1), m, fc and uj 
are the angular oscillation mode, longitudinal wave num- 
ber and angular frequency of an elementary excitation of 
the nanotube atoms valence electrons treated as an elec- 
tron gas. Note that one may set tpo = because of the 
axial symmetry of our model for nanotube. The Fourier 
transform of the induced potential is given by [s^ 

^ind(?') = -[^(''i fij ^) + b^g{r,b,m,k)g'{b,a,m,k)] 

£nt 

+ 63?ff(r,6,m,fc)^|b, (9) 

where g'{r,r' ,m,k) = ■^g{r,r' ,m,k), and the Fourier 
transform of the induced electron number density on the 
nanotube is 



V^$ind(^,t^) = 47r[7'ii(fa,t^)(5(r - a) - crb(fb, cj)(5(r - 6)], 

, , (5) 

where Ta is a position on the nanotube of radius a given 
by = {a, z}, fb is a position on the boundary of the 
dielectric of radius b given by rb = {b, ^p, z}, ni{r^, uj) is 
the induced electron number density (per unit area) on 
the nanotube, and fj\^{f\^,uj) is the polarization charge 
density (per unit area) induced on the boundary of the 
dielectric. 

We may denote the Fourier transform in cylindrical 
coordinates of an arbitrary function A{r, Lp, z, uj) by 

A{r, ip,z,Lu)=J2 [ -^e™'^e^'=^i(r, m, k, u). (6) 

m (271") 

The Green's function in cylindrical coordinates is then 

1 \ ^ f im(ip-ip') ik(z~z') i / _ ; \ 

r= / / 96 ^'e ^ 'q(r,r,m,k), 

\r-P\\ i^nf 

(7) 

where g{r^ r', m, k) is the radial Green's function, 

g{r,r',m,k) = 47r/,n(fcr<)ii'„i(fcr>), 

with r< ~ min{r, r'}, r> — max{ r,r'}, and and 
being the modified Bessel's functions of the first and sec- 
ond kind, respectively. 



l>ext(a) + &.9(a,6,m,fc)5R^|b 

Til . 

^-1 _|_ [^(a, a, m, k) + b^g{b^ a, m, k)g'{h^ a, m, A;)] 

(10) 

Here we have defined the response function of the po- 
larization charge due to the external electric field in the 
radial direction on the outer boundary of the medium, 
with dielectric function by 



47r[£„t + (£nt - e^)kbIm{kb)K'^{kb)] ' 

The response function, x, for the induced electron den- 
sity on the nanotube due to the total electric potential, 
defined by hi — x'^i is given in the 2D hydrodynamic 
model by 

a{k'^ -t- m? I a^) + f3{k^ + m? j a^'f — + 17) 

(12) 

where ng is the equilibrium number density of all four 
carbon valence electrons (tiq ~ 0.428 a.u.), a — nnQ, (3 
= 1/4, and we take the limit 7^0+ [37j . 

We note that the above theory yields the self energy 
(image potential), i?soif(''o)j as a stationary function in 
the proton's moving frame of reference, which only de- 
pends on proton's radial position within an infinitely 
long nanotube. In that respect, it should be possible 
to limit our considerations to nanotubes which are long 
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enough so that we may ignore the edge effects on the 
image potential at the entrance and exit planes and, at 
the same time, short enough that the total energy losses 
of channeled protons may be neglected. Therefore, when 
calculating the dynamic image force on a proton one may 
safely consider its longitudinal velocity component as 
constant equal to its initial speed v, while changes in the 
perpendicular components of the proton velocity u can 
be neglected under channeling conditions. Consequently, 
one may consider the radial position r of a channeled pro- 
ton to evolve adiabatically under the action of an axially 
symmetric force field, F{r) — Frcp{r) + Fimageir), where 
^'rop(?') = -dUrcp{r)/dr and -Flmago(?') = -d-Eseif (r)/dr. 
Because our numerical calculations of proton trajecto- 
ries will be executed under the assumption of a homo- 
geneous, mono-directional beam of protons incident in 
the direction parallel to the nanotube axis, the resulting 
simulation code will be essentially two-dimensional (2D). 
The effects on proton channeling coming from a diver- 
gent beam and non-parallel incidence will be studied by 
a full 3D code in future reports. 

It should be also mentioned that the assumption that 
proton charge remains fixed at Zi — 1, which is used 
in our channeling simulations, needs careful considera- 
tion. Making analogy with the analysis of experiments 
on grazing scattering of protons on an Al surface [46] and 
referring to the data on proton transmission through car- 
bon foils [131, one may expect that the Zi = 1 charge 
state would be the most dominant fraction for protons 
channeled in carbon nanotubes at the speeds in excess of 
V = 3. Nevertheless, simultaneous measurements of the 
angular distributions of transmitted particles in different 
charge states could reveal a wealth of information on both 
the image interactions and the charge transfer processes 
taking place in short nanotubes, in close analogy with the 
experiments done by Winter (Sa . |46| . On the other hand, 
for ion channeling through longer nanoubes, the image 
force itself can be quite strongly affected by the dynam- 
ics of the charge-changing events near the nanotube wall, 
as was shown recently for grazing scattering of protons 
from an Al surface [43]. The effects of charge transfer 
will be included in our future ion channeling simulations. 

The angular distributions of transmitted protons are 
generated using a Monte-Carlo computer simulation 
method. The Cartesian components of the proton im- 
pact parameter, xq and yo: are chosen from a random 
uniform distribution within the nanotube cross-sectional 
area. With the bond length between the nanotube atoms 
being 0.144 nm [4§], we obtain for the radius of the (11, 9) 
nanotube of a = 0.689 nm. Any protons with an im- 
pact point inside the annulus with radii in the interval 
[a — Osc, a], where a^c = [Qtt^/ (128 Z2)] a ao is the screen- 
ing length (with being Bohr's radius), are treated 
as if they were backscattered and are disregarded from 
the simulation. The initial number of protons used is 
3 141 929. 

The Cartesian components of the proton scattering an- 
gle, i.e., of the deflection function, 8x and Qy, are ob- 



tained via the expressions = "x/w and Qy = Uy/v, 
where and Uy are the Cartesian components of the 
final perpendicular velocity vector, u, obtained in our 
simulation of transmitted protons. It has been demon- 
strated that the proton channeling in nanotubes can be 
analyzed successfully via the corresponding mapping of 
the impact parameter plane, the xoyo plane, to the scat- 
tering angle plane, the 0x0y plane. However, as the 
total interaction potential in the case under considera- 
tion is axially symmetric and the incident protons are 
moving parallel to the nanotube axis, the analysis of the 
mapping may be reduced to the analysis of the scatter- 
ing angle, 9 = (9x + ©y)^, as a function of the impact 
parameter, tq = {xq + j/o)^. Therefore, we may take 
j/o = 0, and analyze the deflection function 8x(a;o) only. 
The extrema of this function are the rainbow extrema, 
and the corresponding singularities appearing in the an- 
gular distribution of channeled protons are the rainbow 
singularities 6]. 



III. RESULTS AND DISCUSSION 

In this section we shall first discuss in Figures [1] and 
[2] respectively the image force and the total interaction 
force for a proton channeled in an (11, 9) SWNT in vac- 
uum and encapsulated by Si02, AI2O3 and Ni channels. 
Subsequently, we shall analyze the dynamic polarization 
effects on the angular distributions of protons after chan- 
neling at a speed of ?; = 5 a.u. through an (11, 9) nan- 
otube placed in vacuum and in Si02 for three nanotube 
lengths, L = 0.1, 0.3, and 0.5 /xm in Figures [H |4] and 
m respectively. This analysis will be followed by a com- 
parison of the dynamic polarization effects for different 
proton speeds but with a fixed dwell time, T — L/v = 0.1 
/im/uB — 0.04571 ps (with vb being the Bohr speed), for 
nanotubes in vacuum and in Si02. Next, we shall com- 
pare the effects of different media, namely Si02, AI2O3, 
and Ni, for the case of L = 0.8 /zm, v = 8 a.u., as shown 
in Figures [HI [3 and [HI respectively. Finally, we shall com- 
ment on the zero degree focusing (ZDF) effect for proton 
speeds v — 3 a.u., v = 5 a.u. and v — 8 a.u. for nan- 
otubes of various lengths, placed in vacuum and in Si02, 
as shown in Figure [9| Qualitative analysis of the results 
will also be presented in terms of typical proton trajec- 
tories in the nanotube in vacuum and in Si02, shown in 
Figures [TUj and [ni 

Figure [1] gives the image force of a proton traveling 
paraxially at r = 3 (lower), 7 (middle) and 10 a.u. (up- 
per), with speed v, due to an (11, 9) SWNT in vacuum 
(dashed curves) and encapsulated by three dielectric ma- 
terials (solid curves). The radius of this nanotube is a « 
13.01 a.u. while the radius of the dielectric channel is 
b ^ a + 3.21 a.u. « 16.22 a.u., where we have used car- 
bon's van der Waals radius (« 3.21 a.u.) to approximate 
the separation between the nanotube and the nearby di- 
electric surface. In the case a nickel channel surrounding 
the nanotube, DFT based structural minimization cal- 
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culations performed using SIESTA [4i| for graphene on 
a Ni(lll) surface also yielded a separation of approxi- 
mately 0.17 nm. In Fig. [l]^a) we describe the surrounding 
silicon dioxide (Si02) by a dielectric constant of 3.9 [13], 
in Fig. [2b) we model the dielectric response of anodic 
aluminum oxide (AI2O3) following the method described 
in Ref. [13, and in Fig. [Ijc) we model the dielectric 
response of Ni metal following the method described in 
Ref. [m. The solid and dashed curves in Figs. [1] indicate 
that, for proton speeds below 1; « 3 a.u., the dielectrics 
have little influence on the image force, for proton speeds 
between v ~ 3 a.u. and v ^ 6 a.u. the image force is 
somewhat smaller in the presence of dielectrics, while, 
for speeds above v ~ 6 a.u. the image force is increased 
compared to the vacuum case. These results may be 
explained by the plasmon hybridization in the nanotube- 
dielectric system, as previously described in Ref. 37. In 
brief, protons moving at speeds below w « 3 a.u. do not 
induce plasma oscillations in the nanotube, so that the 
dielectric media are completely screened from the ion and 
do not influence the ion's self energy. Above this speed 
plasma oscillations appear, giving incomplete screening 
of the dielectric media by the nanotube via strong plas- 
mon hybridization. As a result, dielectric media strongly 
affect the proton's self energy and the resulting image 
force at intermediate proton speeds. At high speeds, the 
nanotube becomes increasingly transparent to the pro- 
ton, so that the dynamic-polarization effects tend to be 
dominated by the surrounding dielectric media. This is 
further elaborated in Fig. 1(c), where we also show the 
image force on proton at the same three distances in- 
side Ni channels containing no encapsulated nanotube. 
Specifically, we display the results for a channel of radius 
b — 16.22 a.u. (dotted lines) showing how dynamic po- 
larization of nickel becomes prominent in the combined 
nanotube-channel system only at high speeds. For the 
sake of comparing different channeling systems, we also 
show in Fig. 1(c) the results for image force in a Ni chan- 
nel with its radius equal to that of a (11, 9) nanotube, 
b = a = 13.01 a.u. (thin solid lines). One can see that 
the image forces inside the nanotube in vacuum and in a 
Ni channel of the same radius have generally comparable 
magnitudes, with the case of a Ni channel undergoing 
stronger polarization at high proton speeds, say u > 6 
a.u., and the nanotube providing a stronger image force 
at the intermediate to low proton speeds, say v oc 2-5 
a.u. 

We further analyze in Fig. [2] the total force on proton 
channeling, F = i^iop -I- ^Imagc as a function of the po- 
sition X (in a.u.) across the nanotube radius, for proton 
speeds of (a) u = 3 a.u., (b) u = 5 a.u., and (c) v = 8 a.u. 
Here, the dashed curves denote the case of an (11, 9) 
nanotube in vacuum, while the solid curves of various 
thicknesses denote the effects of the three surrounding 
media, Si02, AI2O3 and Ni, by thin, medium, and thick 
solid curves, respectively. The lower branches of these 
curves show the behavior of the image force i^imagc alone 
close to the nanotube, while the upper branches show 
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FIG. 1: The image force of a proton traveling paraxially at 
three distances from the nanotube axis: r = 3 (lower), 7 (mid- 
dle) and 10 a.u. (upper), as a function of proton speed v in 
a.u., due to an (11, 9) SWNT in vacuum (dashed curves) and 
encapsulated by an (a) Si02 channel, (b) AI2O3 channel and 
(c) Ni channel (solid curves). The radius of the nanotube is 
a = 13.01 a.u. and the radius of the dielectric channel is 6 = 
16.22 a.u. In panel (c) we also show the image force for the 
same three distances inside Ni channels without encapsulated 
nanotube for two channel radii: b = 13.01 a.u. (thin solid 
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FIG. 2: The image force of a proton channeled in an (11, 9) 
SWNT in vacuum (dashed curve) and encapsulated by Si02, 
AI2O3 and Ni channels (thin, medium, and thick solid curves) 
as functions of the proton position x in a.u. across the nan- 
otube radius, for three proton speeds: (a) n = 3 a.u., (b) v = 
5 a.u. and (c) n = 8 a.u. The dotted curve gives the repul- 
sive force, frep, while the lower branches of solid and dashed 
curves give the attractive image force alone, Jlmage, and the 
upper branches give the total force f , due to an (11, 9) SWNT 
in vacuum (dashed curve) and encapsulated by the three chan- 
nels (solid curves). The radius of the nanotube is a = 13.01 



the total force F as a result of adding the bare repulsive 
force i^rep (shown by dotted lines) to the image force. 
Although close to the nanotube the total force is obvi- 
ously dominated by the repulsive force due to the carbon 
atoms, the effects of the dielectric media on the attractive 
interactions due to the image force are greatly affected 
by both the proton speed and the type of the surround- 
ing material. Figure [^a) shows practically no effects of 
the dielectric media for a proton speed of f = 3 a.u., as 
expected from Fig. [TJ At an intermediate proton speed 
of w = 5 a.u, while the curves in Fig. [DJb) for AI2O3, 
Ni and vacuum are very close to each other, especially in 
the inner part of the nanotube, the case of Si02 displays 
a surprisingly weaker image force almost up to the nan- 
otube wall. On the other hand, for w = 8 a.u., the image 
forces shown in Fig. [D^c) have similar values for all three 
dielectric media, especially in the inner part of the nan- 
otube, whereas the image force for the case of nanotube in 
vacuum is substantially weaker at all distances. As com- 
mented in reference to Fig. [1] the nanotube becomes in- 
creasingly transparent at high proton speeds, so that any 
dynamic-polarization effects will be due to the polariza- 
tion of the surrounding dielectrics [37]. This implies that, 
while the image interaction will be heavily suppressed at 
high proton speeds for the nanotube in vacuum, it can 
remain operational in the presence of surrounding media 
and, in fact, will be dominated by their dielectric prop- 
erties as if carbon nanotube was not there. Thus, since 
the image force is crucial for the appearance of rainbows 
and the ZDF in ion angular distributions after channeling 
through short carbon nanotubes, we expect strong effects 
of the surrounding media, especially at higher (but still 
non-relativistic) speeds. 

In Figures we show the results of proton chan- 
neling through the nanotube in vacuum (dashed curves) 
and encapsulated by a Si02 channel (solid curves), by 
both the deflection function QJ^xq) in panels (a), and 
the corresponding angular distributions in panels (b) . In 
Fig. [21 where the proton speed is w = 5 a.u. and the 
nanotube length is L = 0.1 /im, we find a pair of very 
shallow extrema, labeled 1, of the dashed curve in panel 
(a), whereas the solid curve does not exhibit such ex- 
trema. As a consequence, the corresponding angular dis- 
tributions in Fig. [3^b) display, besides massive central 
peaks, also two very small rainbow peaks, also labeled 
1, in the case when the nanotube is in vacuum (dashed 
curve). However, we find no rainbow peaks when the 
nanotube is surrounded by Si02 (solid curve). On the 
other hand, the yield in the central maximum for the 
Si02 case is found to be almost three times larger than 
in the case of a nanotube in vacuum. These findings can 
be explained by examining the results in Fig. [Hb) shown 
by the thin solid curve for the Si02 case in comparison to 
those shown by the dashed curve for a nanotube in vac- 
uum. The rainbow effect is missing in the case of Si02 
because the image force is too weak compared to the vac- 
uum case, and the protons are not pulled far enough from 
their initial direction in such a short nanutobe. On the 
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other hand, the image potential well (not shown here) is 
found to be shallower and broader in the case of Si02 
implying that a larger fraction of the incident protons 
will be very near the initial direction, giving more flux 
of undeflected particles than in the vacuum case. In- 
creasing the nanotube length to L = 0.3 /im, as shown 
in Fig. m amounts to a shift of the extrema 1 in the 
vacuum case (dashed curves), and also to the appear- 
ance of a pair of extrema, labeled by Id, in the solid 
curves due to the image force when the nanotube is in a 
dielectric medium (Si02). A further increase of the nan- 
otube length to L = 0.5 fj,m, as shown in Fig. [5l gives 
rise to multiple rainbows when the carbon nanotube is 
in vacuum. As discussed in Ref. [sl, the rainbows la- 
beled by 1 belong to the first class (rainbow trajectories 
which have experienced one deflection within the image- 
generated potential well), while the rainbows labeled by 
2' and 2" belong to the second class (rainbow trajecto- 
ries which have experienced two deflections within the 
potential well). However, no new rainbow peaks appear 
in the solid curves shown in Fig. [5] when the nanotube 
is in Si02, besides the pair of peaks labeled by Id, which 
are now shifted in comparison to those in Fig. 3) This 
inability of the image force to create multiple rainbows 
in the case of the surrounding Si02 channel can again be 
explained by comparing the thin solid line for Si02 with 
the dashed line for a nanotube in vacuum, shown in Fig. 
Hb). 

We further concentrate on the velocity dependencies in 
the effects of dielectrics on the image force, while elim- 
inating the cumulative effects on ion deflection due to 
the increasing nanotube length. This is done by looking 
into various cases with fixed dwell time, T = L/v = 0.1 
/zm/uB, for nanotubes in vacuum and in Si02. Consid- 
ering low proton speeds, say below u « 3 a.u., we note 
that the results for all three surrounding media (Si02, 
AI2O3, or Ni) should be virtually identical to the results 
for an (11, 9) nanotube in vacuum. This is because pro- 
tons at such speeds do not induce plasma oscillations in 
the carbon nanotube, so that the influence of the dielec- 
tric media is completely screened, as shown in Figs. [T] 
anddja). The results for the combination of parameters 
L = 0.3 /im, V — i a.u. are displayed in Figure 4 of Ref. 
0, where one can see five pairs of the extremal points, 
labeled by 1, 2', 2", 3' and 3", which result from one 
(1), two (2' and 2") and three (3' and 3") deflections of 
the rainbow trajectories within the image-generated po- 
tential well. For the combination of parameters L = 0.5 
/im, V = b a.u., as shown in Fig. [51 one notices a signifl- 
cant depletion of the number of rainbow peaks, down to 
three pairs (1, 2' and 2") for the nanotube in vacuum, 
and only one pair (Id) for surrounding Si02. The further 
reduction of the number of rainbow peaks with increasing 
proton speed is illustrated in Fig. [5] for the combination 
of parameters L — 0.8 /im, v — 8 a.u. Here, we notice 
that the massive central peak is wider than in Fig. [51 
while the rainbow effect has completely disappeared for 
the nanotube in vacuum. This is due to the diminished 
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FIG. 3: The (a) deflection functions and (b) corresponding an- 
gular distributions of protons channeled in an (11, 9) SWNT 
in vacuum (dashed curve) and encapsulated by a Si02 channel 
(solid curve). The proton speed is n = 5 a.u., the nanotube 
length is L = 0.1 /ini, the nanotube radius is a = 13.01 a.u., 
and the dielectric channel radius is b = 16.22 a.u. The angular 
distribution's bin size is 0.0213 mrad. 



image force on protons at such a high speed. On the 
other hand, the one rainbow peak (Id) from Fig. [51 has 
remained in Fig. [6lfor the case of a nanotube surrounded 
by Si02, although this peak is now very small. This per- 
sistence of the rainbow peak for a surrounding dielectric 
can be explained by the "transparency" of nanotubes at 
high proton speeds [131, where the image force is dom- 
inated by the polarization of the surrounding dielectric, 
as shown in Fig. [H^c). 

We also compare the case of proton speed v — % a.u. 
and nanotube length L — 0.8 /tm for the nanotube in 
Si02, AI2O3 and Ni, as shown in Figures [SI [3 and [51 
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FIG. 4: The (a) deflection functions and (b) corresponding an- 
gular distributions of protons channeled in an (11, 9) SWNT 
in vacuum (dashed curve) and encapsulated by a Si02 channel 
(solid curve). The proton speed is w = 5 a.u., the nanotube 
length is L = 0.3 fim, the nanotube radius is a = 13.01 a.u., 
and the dielectric channel radius is 6 = 16.22 a.u. The angular 
distribution's bin size is 0.0213 mrad. 



FIG. 5: The (a) deflection functions and (b) corresponding an- 
gular distributions of protons channeled in an (11, 9) SWNT 
in vacuum (dashed curve) and encapsulated by a Si02 channel 
(solid curve). The proton speed is n = 5 a.u., the nanotube 
length is i = 0.5 /im, the nanotube radius is a = 13.01 a.u. 
and the dielectric channel radius is b = 16.22 a.u. The angular 
distribution's bin size is 0.0213 mrad. 



respectively. One notices that the positions of the rain- 
bow extrema, labeled by Id, hardly change for the solid 
curves shown in Figs. [B][5]for different dielectrics. This 
can be explained by the relatively close proximity of 
the three image- force curves for nanotubes in dielectrics, 
shown in Fig. ^c) for all three dielectrics, at distances 
where extremal points occur in the corresponding deflec- 
tion curves, as shown in Figs. [MH 

As evidenced in Figs. [3][8l the effects of dielectric me- 
dia on the dynamic polarization of carbon nanotubes also 
affect quite strongly the central peaks in the angular dis- 
tributions of channeled ions. We therefore analyze next 



how these effects change the ZDF in the case an (11, 9) 
nanotube in vacuum and in a Si02 channel. We first note 
that our simulations of proton channeling through carbon 
nanotubes without the image force, when proton trajec- 
tories are governed only by the repulsive Doyle- Turner 
potential, did not yield any periodic peaking of ion direc- 
tions parallel to the nanotube. This is explained by the 
short range of the Doyle- Turner potential, which is very 
steep near the nanotube walls, so that protons with dif- 
ferent impact parameters undergo transversal oscillations 
with a wide range of periods. Therefore, there is no single 
frequency of such oscillations that may give rise to the pe- 
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FIG. 6: The (a) deflection functions and (b) corresponding an- 
gular distributions of protons channeled in an (11, 9) SWNT 
in vacuum (dashed curve) and encapsulated by a Si02 channel 
(solid curve). The proton speed is « = 8 a.u., the nanotube 
length is L = 0.8 /im, the nanotube radius is a = 13.01 a.u. 
and the dielectric channel radius is fe = 16.22 a.u. The angular 
distribution's bin size is 0.0213 mrad. 



FIG. 7: The (a) deflection functions and (b) corresponding an- 
gular distributions of protons channeled in an (11, 9) SWNT 
in vacuum (dashed curve) and encapsulated by an AI2O3 
channel (solid curve). The proton speed is v — 8 a.u., the 
nanotube length is L = 0.8 /im, the nanotube radius is a = 
13.01 a.u. and the dielectric channel radius is 6 = 16.22 a.u. 
The angular distribution's bin size is 0.0213 mrad. 



riodic peaking of ion directions parallel to the nanotube 
when their length increases. When the image potential is 
included in our simulations, one finds parabolic regions 
in the total potential near the minima. This allows for 
a broader range of proton impact parameters that would 
give oscillations in the transversal directions at almost 
the same frequency. When the trajectories of such pro- 
tons become almost parallel to the entrance beam, one 
finds the effect of ZDF. In our simulations, we use as 
a criterion for the ZDF that the proton speeds in the 
transverse directions are Vt < 10~^ a.u. 

Figure[5]illustrates the effects of a surrounding medium 



on the ZDF by showing the dependence of proton yield 
(for Vt < 10~^ a.u.) on the dwell time, L/v, for protons 
channeled in an (11, 9) SWNT in vacuum and encapsu- 
lated by a Si02 channel, where the proton speeds are (a) 
V = 3 a.u., (b) u = 5 a.u. and (c) v = 8 a.u. While 
the surrounding dielectric does not effect the ZDF for a 
proton speed of w = 3 a.u., we notice strong effects on 
the ZDF for proton speeds v = 5 a.u. and v — 8 a.u. 
when compared to the case of the nanotube in vacuum. 
Specifically, the period of the ZDF occurrence is seen to 
increase steadily with increasing proton speed in the lat- 
ter case, owing to the weakening image interaction with 
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FIG. 8: The (a) deflection functions and (b) corresponding an- 
gular distributions of protons channeled in an (11, 9) SWNT 
in vacuum (dashed curve) and encapsulated by a Ni channel 
(solid curve). The proton speed is « = 8 a.u., the nanotube 
length is L = 0.8 /im, the nanotube radius is a = 13.01 a.u. 
and the dielectric channel radius is 6 = 16.22 a.u. The angular 
distribution's bin size is 0.0213 mrad. 



the nanotubcs in vacuum, whereas this increase seems to 
saturate at higher proton speeds for nanotubes in a Si02 
channel. Realizing that the ZDF peaks appear more fre- 
quently when the image force is stronger, one can easily 
refer to Figs. [2] (b) and (c) to explain why the first ZDF 
peaks for the nanotubes in Si02 appear in Figs. M (b) 
and (c), respectively, after and before the first peaks for 
the nanotubes in vacuum. 

As regards the yields in the central maxima on the (b) 
panels of Figs. [Uni we have seen for the case oi v = 
5 a.u. and L = 0.3 ^m, shown in Fig. Sfb), that the 
yield when the nanotube is in vacuum is about twice the 
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FIG. 9: Dependence of the zero-degree focusing (ZDF) versus 
dwell time L/v (in units of ^m/vh = 0.4571 ps) for protons 
channeled in an (11, 9) SWNT in vacuum (dashed curve) and 
encapsulated by a Si02 channel (solid curve). The longitu- 
dinal proton speeds are (a) ii = 3 a.u., (b) v — 5 a.u., and 
(c) n = 8 a.u., while the proton yield at ZDF is obtained by 
collecting all transverse proton speeds such that Vt < 10~^ 
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yield when the nanotube is in Si02. This is consistent 
with the fact that, in Fig. [IJb), we are very close to 
the first ZDF peak in vacuum for this proton speed and 
nanotube length. On the other hand, in the case v — 5 
a.u. and L — 0.5 /ini, shown in Fig. [UJb), and in the 
case V = 8 a.u. and L =0.8 /im, shown in Fig. [Hb), we 
notice that the yields of the central maxima when the 
nanotubes are in Si02 are about three times the yields 
for the nanotubes in vacuum. This can now be explained 
by noticing that, for these combinations of proton speed 
and nanotube length, the yields are very close to the 
first ZDF peaks in, respectively. Figs. [5] (b) and (c) for 
nanotubes in Si02. 

The results presented so far may be further elucidated 
by considering the effects of the surrounding medium on 
typical proton trajectories in carbon nanotubes, as dis- 
played in Figs. [To] and [TlJ We note that, in the case 
without image potential, there is only one type of proton 
trajectory, arising from proton oscillations between the 
opposite sides of the nanotube wall. In this case, the an- 
gular distributions of channeled protons through different 
(11, 9) nanotubes are similar for the same dwell times 
L/v, and their widths obey the law v ■ — constant 
for different protons speeds. When the image potential 
is included, there are two characteristic types of proton 
trajectories, arising from proton oscillations between the 
opposite sides of the nanotube wall, and proton oscil- 
lations in the potential minima generated by the image 
potential. 

Figure [TO] shows dependence of the proton deflection 
angle 0x multiplied by the proton speed v on the dwell 
time L/v for an impact parameter of xq = 12 a.u. and 
three different proton speeds (a,) v — 3 a.u., (b) v = 5 a.u. 
and (c) f = 8 a.u. for a nanotube in vacuum and encap- 
sulated by Si02. Since protons with impact parameters 
|a;o| > 11 a.u. have enough transversal energy to make 
oscillations between the opposite sides of the nanotube 
wall, the results shown in Fig. [TO] resemble the motion 
of a particle in a box with rigid walls owing to our use 
of the Doyle- Turner potential. The amplitude of the os- 
cillations incrises slightly due to weakening of the image 
forces with increasing proton speed. More importantly, 
we see very little influence of the dielectric surrounding 
on the type of trajectories shown in Fig. 1101 even at 
higher proton speeds or after multiple oscillations of in- 
side the nanotube. 

Figure [11] shows the dependence of the proton deflec- 
tion angle multiplied by the proton speed v on the 
dwell time L/v for two proton impact parameters, xq = 
3 a.u. and xq = S> a.u, and three protons speed (a) w = 3 
a.u., (b) V — 5 a.u. and (c) v — 8 a.u., for a nanotube in 
vacuum and encapsulated by Si02. Protons with impact 
parameters |a;o| < 11 a.u. do not have enough transver- 
sal energy to bounce between the opposite sides of the 
nanotube wall, but rather undergo oscillations in the 
transverse plane around the minima generated by the 
image potential near x « 9.3 a.u. It is obvious from Fig- 
ure [TT] that the amplitudes of these oscillations decrease 
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FIG. 10: Dependence of the proton deflection angle Ox in 
mrad multiplied by the proton speed v in a.u. on the dwell 
time L/v (in units of //m/ub = 0.4571 ps) for a proton impact 
parameter of xo = 12 a.u. and the three protons speeds (a) v 
= 3 a.u., (b) V = 5 a.u. and (c) u = 8 a.u., due to an (11, 9) 
SWNT in vacuum (dashed curve) and encapsulated by Si02 
channels (solid curve). 
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and their periods increase with increasing proton speed, 
which is easily explained by the weakening of image force 
as shown in Figures [1] and [21 On the other hand, such 
changes in the image interaction did not give any sub- 
stantial effects in Figure [10] because the total potential is 
dominated by the repulsive interaction for proton's im- 
pact parameters |a;o| > 11. a.u., as is obvious from Figure 
[5] However, more important properties seen in Figure [TT] 
are the strong effects of dielectric media on oscillations 
in the proton deflection angle at the speeds w = 5 a.u. 
and V — 8 a.u. Namely, one notices in Figures [TTjfb) and 
(c) that the shapes of oscillations for nanotubes in Si02 
do not change much between those two speeds, whereas 
the periods of oscillations for nanotubes in vacuum have 
much shorter periods for v = 5 a.u., and much longer pe- 
riods for V = 8 a.u. than those in the Si02 cases. These 
notions can be related to the appearance of ZDF peaks 
in Figures [51 (b) and (c) , as well as to the values of forces 
shown in Figures [D (b) and (c) for the same systems. 
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FIG. 11: Dependence of the proton deflection angle Ox in 
mrad multiplied by the proton speed v in a.u. on the dwell 
time Ljv (in units of /im/ub = 0.4571 ps) for proton impact 
parameters of a;o = 3 a.u. and a;o = 8 a.u. and three proton 
speeds (a) u = 3 a.u., (b) u = 5 a.u. and (c) n = 8 a.u., due to 
an (11, 9) SWNT in vacuum (dashed curve) and encapsulated 
by Si02 channels (solid curve). 



IV. CONCLUDING REMARKS 

We have presented the first theoretical investigation of 
the effects of dynamic polarization of the carbon atom's 
valence electrons on the angular distributions of protons 
channeled in an (11, 9) SWNT surrounded by different 
dielectric media. Proton speeds between 3 and 10 a.u., 
corresponding to energies of 0.223 and 2.49 MeV, have 
been chosen, with the nanotube's length varied between 
0.1 and 1 /^m. 

We have confirmed here our earlier findings that in 
short chiral SWNTs, it is the dynamic image interac- 
tion that gives rise to the rainbow effect in the angular 
distributions of protons channeled at the speeds below 
some V K, 8 a.u. @. In the presence of dielectric media, 
this range of the image interaction effects is expanded to 
higher proton speeds because then the SWNTs become 
" transparent" to the polarization effects of the surround- 
ing material. On the other hand, for proton speeds below 
V ~ "i, a.u., the image interaction is almost unaffected by 
the surrounding media because of the efficient screening 
by the nanotube. 

Specifically, we have found that the effects of dielectric 
media are not only of quantitative nature in affecting the 
positions of the rainbow peaks, but rather can give rise to 
qualitative differences compared to the case of nanotubes 
in vacuum, e.g., in yielding different numbers of rainbow 
peaks for the same proton speed and nanotube length. 
For example, we have seen that the presence of dielectric 
media in our simulations both removes (see Figure [2Jb)), 
and introduces (see Figures [S][5fb)) rainbow singularities 
in the angular distributions of protons channeled through 
nanotubes. Moreover, the type of the surrounding ma- 
terial has also been found to affect details of the proton 
angular distributions. 

Going further beyond our previous study Q, we have 
analyzed here the zero-degree focusing (ZDF) effect and 
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found quite substantial differences between the cases of 
nanotubes in vacuum and in dielectric media when it 
comes to the periods of peaks in the proton yield as a 
function of dwell time through the nanotube. These ef- 
fects have been further analyzed by studying typical pro- 
ton trajectories which revealed strong effects of dielectric 
media. Because it may be easy to measure such features 
of ZDF in future experiments on ion channeling through 
carbon nanotubes, we shell devote a separate study to 
the ZDF in chiral SWNTs as particularly simple proto- 
types of quasi-one-dimensional channeling where the im- 
age interaction produces rich structures in the angular 
distributions of channeled ions. 

All our findings indicate that it is important to care- 
fully consider in future simulations and experiments the 
role played by dielectric media in ion channeling through 
carbon nanotubes at the MeV energies. For example, 
the best well ordered carbon nanotubes have been grown 
in porous dielectric media, such as AI2O3, so that in 
any analysis of ion channeling experiments through such 
structures one may not realistically ignore the influence 
of dielectric media. Moreover, in the experiment by Zhu 
et al. , angular distributions of ions channeled through 
the nanotubes grown inside the pores in AI2O3 were com- 
pared with those coming from ion channeling through 
the pores of pristine AI2O3 f?]. It is therefore impor- 
tant to perform comparative simulations of ion chan- 
neling through channels in different materials, with and 
without carbon nanotubes embedded in them, while tak- 
ing proper account of the dynamic polarization effects. 
While a detailed study of this problem is in progress. 



some qualitative predictions can be made here for chan- 
nels inside conducting materials. Based on recent cal- 
culations of the dynamic polarization of such structures 
[30I [sil . [33 , Wf\ , it follows that the image force would be 
expected to play similar role as described in the present 
work, at least judging by its magnitudes for a nanotube 
in vacuum and for an empty channel in metal having the 
same radius as nanotube (compare thin solid lines with 
dashed hues in Fig[T](c)). However, it is questionable 
whether the atomic structure of walls in such channels 
would be smooth enough to support the rainbow effect 
[1^ [l^l . This leaves the ZDF as probably the most ro- 
bust feature associated with the image force, which could 
be probed with comparable degrees of experimental de- 
tail for ion channeling through both pristine channels and 
carbon nanotubes grown in those channels. We finally 
note that such studies may further elucidate dielectric 
properties of carbon nanotubes in the presence of dielec- 
tric media of relevance to nanoelectronic, such as Si02. 
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